Nonlinear magnetic response of the magnetized vacuum to applied electric field 
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We find first nonlinear correction to the field, produced by a static charge at rest in a background 
constant magnetic field. It is quadratic in the charge and purely magnetic. The third-rank polar- 
ization tensor - the nonlinear response function - is written within the local approximation of the 
effective action in an otherwise model- and approximation-independent way within any P-invariant 
nonlinear electrodynamics, QED included. 
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In Maxwell electrodynamics the superposition principle is true, which reads that electromagnetic fields do not 
directly interact between themselves and may be linearly combined independently. This is not the case in nonlinear 
electrodynamics, wherein only small electromagnetic fields are mutually independent. 

A popular example of a nonlinear electrodynamics in the vacuum is provided by the Born-Infeld model [l| , also by 
a noncommutative gauge theory, in this respect considered, e.g., in 0, Q ■ Many issues of nonlinear electrody- 

namics are thoroughly elaborated in [4|. Another, practically most important example is quantum electrodynamics 
(QED). The reason why it is nonlinear is that an electromagnetic field, say a photon, may create virtual electron- 
positron pairs that interact with this field itself and/or with any other, "external" field. This makes a mechanism 
k>" | that lets electromagnetic fields sense each other. 

oo 

\q ■ scattering. When taken off the photon mass shell the corresponding probability amplitude becomes as a matter of 
fact responsible for the leading nonlinear (cubic) correction to the electric Coulomb field [5| 1 that can be conveniently 
written as 



The well-known nonlinear effect of QED, present already in the vacuum without any external field, is light-by-light 
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Here £ — {q/Anr 2 ) is the standard Coulomb field 2 in Heaviside units produced by the point charge q at the distance 



r, while e and m are the electron charge and mass, a = (e 2 /4ir) — 1/137 is the fine-structure constant. It is generally 
known, and also seen in this equation, that in QED the nonlinearity is determined by the ratio of the electromagnetic 
field to Schwinger's characteristic value (m 2 /e) = 4.4 • 10 13 cgse units that makes 1.3-10 V/cm when one measures 
an electric field, and 4.4 • 10 13 G if a magnetic field is concerned. Electromagnetic fields should be comparable in 
strength to these values in order that the interaction between them migt become essential. The nonlinear correction 
in (p} becomes valuable when one is interested to approach a sufficiently small-sized charge sufficiently close. Say, to 
approach the nucleus of a not too heavy atom as close as a few Fm. On the other hand, electric fields, large in the 
Schwinger scale, up to 10 18 — 10 19 V/cm, occur at the surface of strange quark stars Q, depending on whether 
the matter is in the superconducting state @. For such fields the vacuum is unstable, and the Schwinger effect of 
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1 The authors are indebted to M.I. Vysotsky, who attracted their attention to this result. 

2 The linear response to the applied charge q due to the vacuum polarization known as the Uehling-Silber correction to the Coulomb 
potential Q may be also included into £. 
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spontaneous electron-positron pairs by the vacuum becomes already efficient, which requires a special treatment, see 
the book [9(. We do not consider the corresponding complications in the present paper, however. 

In this paper we are dealing with another nonlinear phenomenon, also associated with strong electric field, namely 
the production by it of a magnetic field: this magneto-electric effect becomes possible if an external magnetic field is 
present. 

The linear correction to the Coulomb field of a charge due to the vacuum polarization in a magnetic field was 
studied earlier ficil - fl2| with the result that the hydrogen ground energy level saturates (Tol? 1 as the magnetic field 
grows, and that a string is formed fiof . Some hints were thereby produced for considering 13 1 interquark potential 
in QCD. The nonlinear (purely magnetic) correction to the field of a charge in a magnetic field to be considered now 
for the first time is based on the known fact that in this case not only the photon-by-photon scattering exists, but 
also the photon splitting into two (also two-photon merging into one). The splitting is enhanced by the strength of 
the external magnetic field as compared to the vacuum case above. It was elaborated in theory [lj] and is thought 



of as being efficient in a pulsar magnetosphere with the magnetic fields above 10 12 G [l5j . essentially contributing 
to the electron-positron plasma production and to the radiation pattern of pulsars. Again, the same as above, when 
taken outside the photon mass shell, the corresponding probability amplitudes become responsible for a nonlinear 
induction of time-independent current (and, hence, of the stationary magnetic field) by static charges or, equivalently, 
by static electric fields created by them. The magnetic field produced by a static charge in an external magnetic 
field is even (quadratic in the lowest order of nonlinearity) with respect to its magnitude and linearly diappears with 
the external field - in agreement with the generalized Furry theorem of Ref. [3| that states that the numbers of 
electric and magnetic legs in every diagram should be each even. It also agrees with this theorem in that there are 
no corrections to the static electric field in the lowest (second-power) nonlinear order. Previously, magneto-electric 
effect was considered in 0, Q for classical noncommutative electrodynamics, and within QED as a linear response to 
a static charge by the vacuum filled with external electric and magnetic fields [lj| . 

In the next Section II, for the most general case of a constant and homogeneous external electromagnetic field, we 
outline the derivation of nonlinear Maxwell equations keeping only the first and the second powers of the electromag- 
netic field living above that external field background, and define a notion of a current, nonlinearly induced by a static 
electric field (or by a static charge) . The nonlinear field equations are served by the second- and third-rank polariza- 
tion tensors. In Section III we restrict the external background to the magnetic-like field, i.e. the one that is purely 
magnetic in a class of special Lorentz frames. Then the involved polarization tensors are given in small-4-momentum 
limit, called also infrared or local approximation, in terms of the derivatives of the effective Lagrange density over 
the background field invariants, bearing in mind that in the local approximation this density does not depend upon 
space-time derivatives of the background field strength. In Section IV we are working in a special frame, where the 
background field is purely magnetic and the static charge is at rest. We calculate the nonlinearly-induced current and 
its magnetic field as expressed through the static electric field produced by the charge. The limiting cases of very 
large and very small background magnetic field are discussed within QED referring to the one-loop Euler-Heisenberg 
effective Lagrangian. In Conclusions the results are resumed, and numerical estimates of the domains of their appli- 
cability are given. Detailed calculations of the second- and third-rank derivatives of the effective action used in the 
work are presented in Appendix within the necessary local approximation. 



II. NONLINEAR ELECTROMAGNETIC FIELD EQUATIONS OVER A CONSTANT FIELD 

BACKGROUND 



In QED and in any other U(l)-gauge-invariant nonlinear electrodynamics, the field equations, when written up to 
terms, quadratic in the small electromagnetic field potential a v {x) 1 have the form 



r) n-d"d v ] a v {x) + / dVlV(x,x , K(.T , ) + 



+ - / d A x'd 4 x"Tl pucT (x 1 x'x")a l/ {x')a° {x") = j p (x), 



(2) 



where j p (x) is a (small) source of the field, Greek indices span the 4-dimensional Minkowski space taking the values 
1,2,3,0, the metric tensor is 77 „ = diag (1, 1, 1, —1), and □ = V 2 — dfi. The second- and the third-rank polarization 
tensors, H pL , and n pl/<T , here are, in the presence of an external field potential A@ (x) = A^ xt (x), defined as 



fl/iT (X, X ) 
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5Av(x)SA T (x') 



(3) 
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8A»(x)5A T (x')5A°(x") 
in terms of the effective action 



(4) 



r = J £(z)d 4 z, (5) 



the generating functional of all-rank polarization tensors - the vertex functions - known in QED as the Legendre 
transform of the generating functional of the Green functions [13] ■ The parameter of the power expansion, to which 
Eq.([2]) provides two lowest terms, depends on a field scale of a definite dynamical theory. We shall discuss this issue 
in Section IV below for QED. 

We did not write the zero-power term (a u (x))°, an external macroscopic current, in equation ([2]), because we 
assumed that the external field had been subjected to the sourceless field equation 



SS 



SAP (y) 



= 0, (6) 



where 

S = jL{z)d 4 z : L(z) = -$(z) + £(z) (7) 

are the total action and the total Lagrangian, respectively. Here — $(z) =(1/4) F^F V ^ is the (free) Maxwell La- 
grangian, F a g{z) — d a Ap(z) — A a (z) is the field-strength tensor. In what follows we shall only deal with external 
fields F a g = d a A° a * t —d l3 A'£ ct , which are independent of the 4-coordinate z^, and with the case where the effective 
Lagrangian £(z) may depend on z M only through the field tensor F a g(z) and its space-time derivatives, and not 
explicitly. The latter property is fulfilled in QED and will be also assumed for other theories subject to our consider- 
ation. Under this assumption the constant field does satisfy the exact sourceless nonlinear field equation ([5]). To see 
this, we fulfill the variational derivative in it 



SS 



5A B (x) 



= 2^ 



SS 



dz 



-^-<5 4(n) (z - z)d 4 z, 



SS 



where (n) marks the derivative with respect to any space-time component. Once the variational derivative SF in)^ ^ 

when restricted onto the coordinate- independent fields F^ u (z) — F lllll cannot depend on z, the integration by parts 
turns this integral to zero. 

The above presentation explains why Eq. ^ is the field equation for small electromagnetic perturbations a°(x) = 
A@ (x) — A^ xt (x) over the external field of a constant field strength, caused by a small external current j p (x) and 
taken to the lowest-power nonlinearity. 

Polarization tensors of every rank II /JT ... CT (a;, x', ...x") satisfy the continuity relations with respect to every argument 
and every index (the transversality property) 

^ r U^.. T ... (T (x,...x',...x")=0, (8) 

necessary to provide invariance of every term in the expansion of T in powers of the field a" under the gauge 
transformation of it. Note that this is the primary property of T as a functional given on field strengths and their 
space-time derivatives only. 

In our case of the external field with space- and time-independent strength the translational invariance holds true, 
which fact makes the all-rank polarization tensors depending on their coordinate differences. 
With the definition of the photon propagator D^(x,x') 

D-}(x - x') = [n^D - <5 (4 V -x) + U^(x - x') (9) 

the nonlinear field equations © take the form of (the set of) integral equations 

a\x) = [ d 4 yD x ?(x - y)j p {y) + [ d 4 yD^(x - y)j*(y), (10) 
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= ~ g J d 4 yd 4 uU^ a (x -u,y- u)a T (y)a°(u), (11) 

where we have introduced the notation j^(x) for what we shall be calling "nonlinearly induced current". 

Before proceeding, the following explanation seems to be in order. Within the present approach the electromagnetic 
field a x (x) is not quantized, this is not needed unless we leave the electromagnetic sector. The nonlinear equations 
written in this section are classical and will be treated classically below in understanding that the effective action is 
known. In QED the latter is the final product of quantum theory, obtained by continual integration over fermions [l7| . 
The effective Lagrangian and all-rank polarization tensors involved are subject to approximate quantum calculations 
and, hence, are functions containing the Plank constant, electron mass and charge. Available is the effective action 
in the local limit referred to in the next section, which is known as the Euler-Heisenberg action when it is calculated 
within the approximation of one electron-positron loop, see @, and as Ritus action when it is calculated with two- 
loop accuracy [18| . The second-rank polarization tensor ([3]) was calculated in the one- loop approximation when the 
external background is formed by a constant and ho mog eneous electromagnetic field of the most general form (when 
the both its invariants 3" and 25 are nonvanishing) in [19j . One-loop diagrams with three photon legs corresponding to 
the third-rank tensor (jj) were calculated both on and off the photon mass shell for QED with external magnetic-like 
(£ > 0, © = 0) and crossed (# = 25 = 0) fields in [14|, and for charge- asymmetric electron-positron plasma without 
external field using the temperature Green function techniques, in [2Cj. The calculations of Stoneham in [14] might 
become a basis for extending the results of the next sections beyound the local approximation used in getting them, 
but they are overcomplicated and not well-structured, so we leave this extension for future. In the next sections we 
stick to the general form of the effective Lagrangian and refer to its specific Euler-Heisenberg form only at the very 
last steps for getting numerical estimates. 



III. LOCAL LIMIT 



From now on we shall restrict ourselves only to slowly varying fields a x (x) and, correspondingly, to consideration 
of the sources j p (y) that give rise to such fields via equations (Tl0|) . (fTT|) . To this end we may take the effective action 
in the local limit, where the field derivatives are disregarded from this functional. This is equivalent to going to the 
infrared asymptotic limit in the second- and third-rank polarization operators, i.e. to keeping, respectively, only the 
second and the third powers of the 4-momentum k p in their Fourier transforms. Aiming at the local limit, we may 
admit that the effective Lagrangian £ depends only on (the relativistic and gauge invariant combinations of) the fields 
F p(7 . Moreover, as long as constant fields are concerned all such combinations may be expressed as functions of the 
two field invariants $ — \F p<J F f " T and 25 = jF pa F pcr , where the dual field tensor is defined as F pa — ^e pr7 \ K F XK 7 with 
the completely antisymmetric unit tensor defined in such a way that £1230 = 1- Then the variational derivatives in ([3]) 
and (j4]) can be calculated in terms of derivatives of £ (S", 25) with respect to the field invariants reduced to the space- 
and time-independent external field. Henceforth we shall be interested in the special case, where the external field is 
a constant purely magnetic field in a certain class of reference frames, called special below. Since in other Lorentz 
frames the electric field is also present we refer to this case as magnetic-like. The invariant conditions that specialize 
the magnetic-like case are 3" > 0, 25 = 0. Once the invariant © is a pseudoscalar, the Lagrangian of a P-invariant 
theory, QED included, may contain it only in an even power. Hence all the odd derivatives of £ (3^, 25) with respect 
to it disappear after being reduced to the external magnetic-like field: 



<925 



d<8d$ 2 



F=J r ,&=0 



F=J r ,0=O 



d 2 £(3,6) 



<925<93 
d 3 £(3, 0) 



d<5 3 



Q 



F=F,®=0 



(12) 



We calculate (|3|) and (|4]) in Appendix using the rule 



SF aP (z) 



9 9 \ s-i, 



6A>*(x) 

(understood as integrated over z with any function of z) by repeatedly applying the relation 



(13) 



ST 



ftC(ff(z),<8(2)) 0£(g(z),0(z)) - - 

- r a P .\ z ) T „ —, r -Tafj.\Z) 



<925(z) 



-^-5\x - z)d i z 
oz a 



(14) 
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and reducing the final results onto the external field. Then, taking eqs. (|T2j) into account and using the notations 



dc(gr,o)) 



-333 : 



d£ 
d 3 £(S,0) 



> £33 



F=J r 



&z 2 



£®e 



3 2 £(S, <») 



d® 2 



F=J r .&=0 



d£ 3 



_ d d 2 £($,<8) 
F= r> ~ d£ d& 



F=.F,©=0 



all with Ffa, = Tfui substituted (hence, from now on, £ = \T p(y T pa > and © = \T pa f p(J = 0), we get for the 
second- 

n^Or - y) = £ 5 (J?^ - v n) 5 4 (x - y)- 
and for the third-rank tensor in the infrared limit 



(15) 



tf^A x -y,x-u) = 

= -0^a a p 1 ^ r ((^ r S 4 (y-x)) (-^S^x-u) 



dx a \ \ dx 



(16) 



where 



(17) 



(The reader may consult Appendix for detailed calculations.) This tensor turns to zero, when there is no external 
field, T — 0, in agreement with the Furry theorem. The two transversality conditions © for (fT5")) are provided by 
that the matrix in the brackets is antisymmetric under each permutation \l -4=> a and r -<=>■ /3, while the first term 
in (fT5l) is transverse explicitly. The three transversality conditions ((HJ for (fT6|) are provided by that matrix (fl"7|) is 
antisymmetric under each permutation /i a, r <^=> /3 and er 7. Thanks to the two latter antisymmertries, 

by using (|16p . (|17[) in (|11[) we obtain for the nonlinearly induced current the expression 



3l\x) = -0^ T _ (f^ r «) , ( i 8) 

that includes only the field intensity tensors / /3r = -^ja T {x) — -^-a^(x). Therefore, the nonlinearly induced current 

is gauge- invariant: it depends only on field intensities and, besides, it is conserved, gf-j'JJ 1 ^) — due to the first 
antisymmetry /1 •<=>■ a. 

We have to approach the nonlinear set (p~0|) . (|TT|) by looking for its solution in a power series in the field a x (x). Within 
the first iteration, to which we shall as a matter of fact confine ourselves, we substitute the linear approximation to 
the solution of equation (jTUJ) 

a}?{x) = J d 4 x'D vp (x - x')f{x') (19) 

for a(x) into (fTTj) . In other words, we should use the electromagnetic field fp T = fj™ = ^^a^ n (x) — ■^p F a 1 p n (x) 
linearly produced by the source j^ix) in the expression for the nonlinearly induced current ()18j) . 
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IV. MAGNETIC FIELD OF A STATIC CHARGE AT REST IN EXTERNAL MAGNETIC FIELD 



We are in a position to start studying the nonlinear effect of production of a magnetic field by a static charge at rest 
in a constant and homogeneous external magnetic field in a special frame. The linear effect of the external magnetic 
field on the electrostatic field of a charge was studied earlier (beyond the infrared approximation) in [Tol - [T2j . 

In this frame the external magnetic field is defined as Bi — (l/2)cijkJ r jk = Fm, B — |B|, while the external electric 
field disappears Ei = J-Qi = 0. The roman indices span the 3D subspace in this reference frame, e^-fc is the fully 
antisymmetric tensor, e 123 = 1. 

Consider now a static charge given in that frame by the 4-current j^{x) — j^(x)5^Q. In the linear approximation 
(1191) , naturally, only an electrostatic field is generated in that frame. Hence, the components with a = f3 = 7 = 0, r, 
a 7^ do not contribute to (|18p. so we need only the components 



O 



iOOjmn 



= -TV 



'66 



CijkBt [ 5 mn £$$ — B n B m 2$@@] + [B m €ji n + B n eji m ] £@q 



(20) 



in (|18p (and the ones obtained from (|2U|) by permutations between the second and the fifth, and between the third and 
the sixth indices), while Oootymn = according to (IT7|) . Therefore jo (x) = 0, i.e. there is no nonlinear (quadratic) 
correction to the static charge within the current quadratic approximation: the induced current (|18[) is purely spacial: 



3 

^iOOjmn \JmOJnO, 



(V X B) i % 5 5 2 - % 06 {B£f - £ ee (V X £\ (BE) , 



(21) 



where £ n = £ n (x) = f { 



lin 

On Ox 



d r .lin 



Q m (x) is the time-independent electric field, linearly produced following Eq. p^)) . and 
the differential operator V acts on everything to the right of it. The magnetic field strength h(x) generated by this 
current according to the Maxwell equation 'Vxh(x) = j n (x) is 

hi(x) =t) l {x)+W i n, (22) 

where 

fo(aO =y [& ss (£(x)) 2 - Sbee (B£(x)) 2 ] - fj(x)£ee (B£{x)) , (23) 
because Vx V£! = 0, and the scalar function O should be subjected to the Poisson equation 

to make the magnetic field h(x) obey the other Maxwell equation 'Vh(x) =0. Hence, the magnetic field is the 
transverse part of (f2"3"| : 



f)i(y) , 3 



\x - y\ 



d 6 y- 



(24) 



Note that the substitution of the field of a point-like charge into (I24|) through (|23|) would cause the divergency of the 
integral in (|2"4")l near y = 0: the present approach fails near the point charge, since it is not applicable to its strongly 
inhomogeneous field. Dealing with the point charge would require going beyond the infrared approximation followed 
to in the present work. Nevertheless, Eq. (|2"41 is sound as applied to extended charges. 

Eq. (|2"4")l would coincide with the magnetic induction b(x) = V X a nl (x) if the linear vacuum magnetization effect 
might be neglected, i.e. if the nonlinear correction to the field in (flOl) 



a nl (x) 



d 4 yD\(x-y)jf(y) 



(25) 



might be taken without the contribution of the linear response function (|16p Yl^^x — x) in the photon propagator ([9J. 
Taking this contribution into account results in more complicated integrals. The situation remains simple, however, 
when we may disregard the anisotropy of the linear magnetic response. The inverse magnetic permeability tensor 
inherent in the second-rank polarization tensor (|15p is, in the special frame, the constant tensor [2l| 
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My = (1 - £$) $ij ~ £$$BiBj, 

whose two 3 eigenvalues fij 1 = 1 — £5, and jUj7 = 1 — £ 5 — 25"£ ss are responsible for magnetizations linearly caused by 
certain conserved constant straight-linear currents flowing along the external magnetic field, and across it, respectively 
(see Appendix in [Hj]). In QED, the values £g and 2$£$g are of the order of the fine structure constant a = 1/137, 
but depend on the field B. When B is very large, B ^> m 2 /e, these quantities, as found from the Euler-Heisenberg 
one-loop effective Lagrangian, behave as, see, e.g., [23| 

a , eB a 
%«— In— , 23%5«— • 

6TT Tfl 6TT 

So, when ^| 3> 2.7, the contribution of 2$£g$ may be neglected as compared to £5, and the linear magnetization 
becomes isotropic, /ij 1 = M^ 1 - Therefore, in this limit, we finally have for the nonlinear magnetic induction 

b(x) = (l-£y)~ 1 h(x). (26) 

The electric field £ = — Voq (x) to be substituted in l|2ip and l|23p is the one that is linearly produced via Eq. 
(|19P by a static charge distribution within the same infrared approximation. To determine it, note that in (|19p 
only the propagator component Doo participates, that, in the Fourier representation, is Dqo = {k 2 — Jfz) - 1 , with >ci 
being one (out of three) eigenvalues of the second-rank polarization tensor ([3]) taken in the static limit feg = in the 
special reference frame. Once the polarization tensor is considered in its infra-red limit (|15p , this quantity is [22I [23( ] 
K2 = k 2 £$ — k 2 2$£®&. Here k 2 = k\ + fcjj, and fcj_, k\\ are the momentum components of the small electromagnetic 

field across and along B , respectively. Now, the calculation of (1T91) for the point-like charge jo{x) — q8 3 (x) results, 
with the use of this propagator, in the anisotropic Coulomb law 

a^(x) = ± 1 (27) 



where e± = 1 — £g and e y = 1 — £5 + 23 r -C<»© are eigenvalues of the dielectric tensor [2l| = (1 — £5) Sij + 
£&@BiBj, responsible for polarizations, linearly caused by homogeneously charged planes parallel and orthogonal to 
B, respectively, x± and xn are the coordinate components across and along B. For large magnetic field one gets 

the linearly growing asymptote from the Euler-Heisenberg Lagrangian 2g , £ a j0 w — jr. This means that if ^4> — 

the dielectric component dominates over e±, i.e. the electrization becomes highly anisotropic, in contrast to the 
magnetization. In this asymptotic region Eq. (|27fl becomes (if we disregard the polarization in e± by setting e± = 1) 
the large-distance behavior of the potential of a point charge in a strong magnetic field calculated in the linear 
approximation in [ToL [ll[ beyond the infrared approximation of the polarization tensor. Note that Eq. (|27p . as well 
as its high-field limit, is only valid far from the charge. In that domain, however, it also fits any charge, with the total 
value q, distributed over a finite region. 



V. SOME NUMERICAL ESTIMATES 



To analyze the large magnetic field limits of the induced current (|2ip , of the resulting magnetic field (f24| and of its 
induction (f2l))) one should also bear in mind the asymptotic behavior £3050 = — 3rr ^i B i ■ Then it follows from the large 
external magnetic field asymptotic behavior, , of the other derivatives of the Euler-Heisenberg Lagrangian 

involved in (l23l) that were listed above that in this limit 
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3 The constant background magnetic-like field makes an uniaxial medium in any of the special frames |2ll . 
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The minus sign in the first line undicates that the induced magnetic field diminishes the external field in the large 
external field regime. 

We may apply the results (|23|) . (|24|) to small external magnetic field (eB/m 2 ) « 1, as well. With the Euler- 
Hciscnbcrg Lagrangian density, one has in this regime: 



4a / e y (a ( ay a / e y 

457T \m A ) Abir \m z I 6l5n \m l I 

The third coefficient £500 does not contribute in the leading order in (eB/m 2 ) << 1 to the estimates 



J " 45tt 



- 7 



e£ii 



7a /e£ ± \ f e£\\\ eB 

U~-B-rr -F 2 ' — <<L 
4o7r \ m z / \ m J m z 

In this approximation we may set e± = e (| = /ij 1 = /i^ 1 = 1. Therefore, h = b, and for the electric field of a charge 

outside of it, one may use here the standard Coulomb law £ = (q/4ir) x/\x\ 3 instead of ([27]) . 

Note that a/457r = 5 • 10~ 5 . So, for the electric field value close to Schwinger's 1.3T0 16 V/cm, the nonlinearly 
produced magnetic field makes up to 3- 1CP 4 of the external magnetic field, which must be kept below Schwinger's 
4.4-10 13 G in this case. 



VI. CONCLUSION 



In this paper we have found an expression for the magnetic field h(x) produced by a static charge q placed into 
an external magnetic field B, Eqs. ([2"3")) . (|24[) . It is shown that in QED this nonlinear magneto-electric effect, not 
considered before, occurs already in the simplest approximation, where the effective Lagrangian £ is taken in its local 
limit, and only second power of the charge q and/or of its electric field £(x) are kept. As for the background magnetic 
field B, to reveal the effect suffice it to take it into account in the linear approximation, ~ B, although magnetic field 
B of arbitrary magnitude is included in our result, as well. The final formulas depend on the first three derivatives 
of the effective action £ with respect to the external field invariants, which complies with the fact that, minimally, 
diagrams with three photon legs are responsible for the effect in the given approximation. 

The results are model- independent and relate not only to QED, but also to any nonlinear electrodynamics provided 
the standard postulates of U(l)-gauge-, Lorentz-, translation-, C-, P-, T- invariances are respected. When applying 
them to QED we take the Euler-Heisenberg Lagrangian for £ to estimate the regimes of weak and strong B. In 
QED all electromagnetic fields appear in ratios to the Schwinger characteristic value m 2 /e of 4.4- 10 13 cgse units. The 
nonlinear magneto-electric effect we are reporting on is efficient, if the electric field of a charge is comparable, but 
still smaller than m 2 /e. Such fields take place near atomic nuclei and at the surface of a strange quark star. Strange 
quark stars can be strongly magnetized, besides [24| ■ When the Schwinger value is exceeded by the electric field, the 
nonlinearity can no longer be treated via the power expansion ^ , and also electron-positron pair creation from the 
vacuum must be taken into account. 
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Appendix 



The second variational derivative of the local effective action is 
5 2 T 



5A»(x)5A T {y) 
, d£($(z),<S(z)) 



d 4 z 



f 3£(g(z),g(z)) , , 

1 ^77i VlnrVaP ~ V M pV aT ) + 



d<8(z) 

2 £(g(z),g(z)) - ~ 
^2 -F«^(2)-F>r(2) H H7^7I\\2 F all (z)F l 3 T {z)+ 



dim? 

9 2 £(S(z),<3(z)) 



(It was convenient to exploit the relation 



<9(©(z)) 2 

F a ^z)Fp r (z) + F atl (z)F 0T (z)] I (^S\x - z)) (£-_fi*(y - z) 



SF a/3 (z) d 4 



(28) 



following from (fT"3|)). 

After reduced to the constant and homogeneous magnetic-like external field this gives the second-rank polarization 
tensor (Tl"5j) . The latter can be also used when the external field is "crossed", i.e. when & = 0, $ = 0, but F pa 7^ 
(electric and magnetic field vectors are mutually orthogonal and equal in length in every Lorentz frame). 

The next derivative of (|2"51) is 



<5 3 r 



8At*(x)5A T (y)5A' T (u) 

ra 2 £(?(z),©(z)) 



9 2 £(3W,S(^)) 



F 1<T (z) + 



9 2 £(S(z),©(z)) - 



dS 2 (z) ' d®(z)d$(z) 



F ia {z) 



d3{z)d<8(z 
d 3 £(3(z),0(z)) 



F ia {z) + 



d 2 £(5(z),©(z)) ~ 



9(©(z)) 2 



F ia {z) 



<9z- 



<5 (u — z)+ 



d*z{3{z)M*)) 
d(3(z)f 

+ F afl (z)Fp T (z) 



F a ,{z)Fp T {z)F ia {z) + ^f^lf^ F^(z)F 0T (z)F^(z) 



F a p(z) T) 



d 
dz~P 



d<&{z)d{Z{z)Y 
d 



dz 1 



d_ 

dz c 



d 



^-5\u-z)- 

OZj 

5 4 (u — z)+ 



a 3 £(g(z),0(z)) d 3 £(ff(z),g(z)) ~ . 



dz(z)d(<s(z)y 



d${z)d<S{z) 



_d 

dz 



d{®{z)f 

■^-5 i (u - z) + 

oz 1 



dz. 



-<5 4 (u — z) + 



g 2 £(ff(z),0(z)) 

a(©(^)) 2 

d 2 £(3(z),©(z)) rr , , . a - , ./ a 3\ - , ./ 9 9 

-[F ali {z)ep Tla — + F [3t {z) [ v ^—- V a *Q^ ) + F a p(z) IVt*q^ ~ Vpaj^ 



dz^ 



'dz" 







+ Ff; T (z)e aflla —}6 4 (u -z) + 



9 f i 
dz- 



F a p{z)Fp T {z) + F a ^{z)Fp T {z) 
6 4 (u-z)} 



d 3 Z(d(z),<5(z)) F ^ (z) | g 3 £(^),®(z)) -^ ( ^ 



a^ 2 (z)50(z) 



9^(z)50 2 (z) 



' ? 5 4 (x-*) > ) (-?-S*(v-z) 



<9z Q 



(29) 



Each z-derivative here and in the next relation is meant to apply only to a single (5-function, the closest on the right 
of it. 

Eq. (|29l) can be used for calculating higher variational derivatives, when needed. 
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After we set 25 =0 and use (TT2|) we get for the P-even theory 



<5 3 r 



8A^(x)8A T (y)8A a {u) 



d 2 z{$(z)Mz)) 



F ia {z) 



+ e, 



d{<5{z)) 2 



d($(z)) 3 

a 2 £(g(z) ; @(z)) 

dim) 2 

+ F af +(z)Fp T (z) 



F a ^z)Fp T {z)F ia {z) 



dS 2 (z) 

8 A (u — z) + 
d 



9 X 4< 

dz~, 



+ 



F a ^{z) rj 



d 



d 



d 3 £(Z(z),<8(z)) 

2 -FycryZ) 



dz(z)d(<&(z)y 



'dz c 



'dz» ) 



8 A {u~z) + 



d 2 z{${z),e{z)) 



F afl (z)ep T ~ ( a + Fp T (z)e 



+ 



d(&(z)) 2 
F a „(z)Fp T (z) + F ail (z)F Pr (z) 



d 3 £(S(2),0(z))~ ' 



d$(z)d<5 2 (z 



It is understood that © is set equal to zero also in the right-hand sides of this and the next equations. 
We integrate by parts with the delta- function S A (x — z) to get 

S 3 T 



0=0 



8A^(x)8A T (y)8A' T {u) 

J d 4 Z(5 4 (x - Z)^~{ (v^rVaP ~ VtfVcrr) 

"5 3 £(^(z),©(z)) 



d 2 £(3(z),<5(z)) 



d(3(z)f 
a 2 £(5(z),«(z)) 



-F a ^{z)Fp T {z)F ia {z) 



d$ 2 (z) 
+ F af ,(z)Fp T (z) 



F ia {z) 

d 3 £($(z),®(z) ) 
d5(z)d{&(z)) 



2 F 1<T (z) 



d(<8(z)) 2 
F ail {z)F 0T (z) + F ail (z)Fp r (z) 



d 2 m{z)Mz))t 



d{<S{z)f 



F ia {z) 



d 3 £($(z),<S(z)) - 

r 1 a\Z) 



d${z)d<S 2 {z) 



d 



dz~, 



-—8 (u-z)\ -—8 (y- z 







d A z8*(x - z) 



8 



r dzJ 



d d 2 Z($(z),&(z)) 

dz a d{z{z)) 2 



'dz c 



'dz» 



d 

8 4 {u- z)—8 A {y - z). 

dzp 



Here and in the next relation the left-most derivative d/dz a acts on all functions of z to the right of it, whereas each 
other z-derivative acts only on the first 8- function placed to the right of it. 

The third-rank polarization tensor (|16[) is obtained from this expression by reducing onto z-independent fields 
F(z)=T: 
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—d 

U ]^a( x -y,X-U) = —{ (v^Vaf) ~ VrfVar) 



Q$2 T l° 



r a 2 £(^,0) 



d& 2 



d$ 3 
d 2 £{$, « 



-^—6 4 (u - x) + 

ax 1 



dx~, 



d 



9\ r / d d \ 



S 4 (u-x)+ 



d3d<5 



T 

2 J 7cr 



dx> 



-5 4 (u — x) + 



d® 2 



_d_ 



a 3 £(g,6) - 



2 J 7 cr 



(5 4 (w-a;)+ 

9^ 54( "- a;)} ^ 4(2/ - a;) - 



It is meant that © = \j- f " J J- pa = 0, and $ — jJ r prT J rf " 7 here. This formula is also applicable to the crossed field 
= 0,^ = 0,^^0. 
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